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1. Introduction 



The perturbation theory in QED used the exact solution with taking into 
. account the special form of interaction is constructed. The mean electro- 

magnetic field of charged particle is calculated. The possibility of elimina- 
tion the problem with ultraviolet as well as infrared divergences is shown. 
The electromagnetic energy of the particle turns out to be regular and 
q ■ small 

Pk 
^— > 
a. 

In the conventional Dyson- Feinman approach in quantum electrodynamics 
all quantities are calculated within the perturbation theory with using the 
zeroth order photon and electron propagators (see, e.g., the famous mono- 
graphs ([1],[2],[3],[4],[5]). Fundamental difficulty of the approach is a diver- 
gence of integrals representing radiative corrections to the electron mass, 
electron charge or power of interaction between electrons and photons. A 
renormalization procedure allows us to eliminate the divergences with the 
help of any form of subtracting the infinitely large terms in the perturbation 
theory. Some famous physicists (including Feinman [6], Landau and Lifshitz, 
Gell-Mann [5]) were estimating this procedure as incorrect. There is a num- 
ber of attempts to solve the problem with divergences in the framework of 
conventional theory. In the last few years some works was published where 
the authors find the possibility to avoid the divergence problem as in QED 
as well as in the other field theories with the help of "clothing" procedure 
firstly proposed in [7] (see, e.g., [8], [9], [10], [11], [12] and references therein) . 
In this connection one should note, that the "clothing" is fulfilled within the 
conventional perturbation theory. On the other hand, it is interesting to 
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know, whether this difficulty is the property only of the perturbation the- 
ory or of the quantum electrodynamics itself? For answer this question one 
should construct any other form of the perturbation theory, quantitative dif- 
ferent from the Dyson-Feinman approach. We would like to show that the 
appropriate form of the theory indeed can be constructed where as infrared 
as well as ultraviolet divergences turn out to be absent. The general scheme 
of the theory was described in [13]. The relativistic units {% = 1, c = 1) are 
used throughout the paper. The normalization volume is set equal to unity. 
We use the common-famous notations of four-dimensional relativistic theory 
of fields. The Greek indices run over 0,1,2,3; the Latin one - 1,2,3; g a p is the 
metric of the pseudo-Euclidean space-time with the signature (1,-1,-1,-1). 

2. Basic Definitions 

Consider the interaction of a free electron with the field of its own radia- 
tion. Consider electron as quantum of the spinor field with the wave operator 
in usual form 

^ = E [uxpe- ipx a X p + vxpe ipx d{p] , (1) 

X,p 

where p - is the momentum of the electron, A = ±1/2 denotes the projection 
of the electron spin on the the specific axe of quantization (down we set this 
one is z-axe). Bispinors taken part in (1) are defined, e.g., in [5], 



UXp 



y/ep +mw\ 



yflep \^ep-m (npa)w X/ 



VXp = U-\-p 



V^p-m {npB)w' x 



Here o - is the vector of Pauli matrices, p° = ep = (j9 2 + m 2 ) 1 / 2 - is an energy 
of the electron obeying the momentum p. The unit spinors w\,w' x describe 
the states with the spin projections equal to A . In the equation (1) d X p - 
is the electron destruction operator in a state with the defined polarization, 
momentum and the energy ep, d\p - is the positron creation operator for 

the same state. All the operators d\p, a X p, d\p, d\p obey the standard Fermi 
commutation relations. 

In the current density operator 
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where :...: means the normal ordering, e - is the electron charge, ip = -^7° 
is the usual Dirac conjugate operator. In the following we need the Fourier 
representation of the current density operator which in our case has the form 

\,a,p 

a A,p / LJ (T-p-q u \p Uj a-p-q e ' / a a-p+q^\p^a-p+q e 

A? e" 4 ^-^] . (2) 

Note, as a rule, the current density operator don't conserves the spin of the 
particle. The conservation of the spin is possible only at q = . 

Now we need to construct the " approximate" current operators which are 
to be commute if related to different times. For the one electron problem one 
of the appropriate expression is given by the formula 

j^)(t) = rmp q , (3) 

The 4- vector / M (<f, t) is defined below and 

\,k 

is the Fourier component of the time-independent electron density operator. 
Note, the operators (3) include only diagonal terms in polarization indices 
and if even refer to different times, obey the commutation relations 

[Jfe(*).5"f(0)(O]- = 0. (4) 

Down we consider (3) as the zeroth order approximation to (2). To define 
the vector / M ((f, t) we impose the condition: the current mean value of (3) 
is to coincide with the mean value of (2), namely, 

(t\jl(t)\t) = nq,t)(t\pt\t). (5) 

This condition ensures the been chosen f^(q,t) gives the sufficiently rapid 
convergence of series in the modified perturbation theory. Obviously, the 
equation (5) may be solved only approximately because the exact vector of 
state is unknown. 

3. Solution to the Ground State Problem 
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Denote the deviation of the current from its " zeroth" value as 

A#(t) =#(«)- $„)(*) 

and represent the electromagnetic interaction as a sum of two parts, Hi n t(t) = 

where 

H\%) = e J jf 0) (r, t)A,(r, t)dV; (6) 

H$(t) = e J Af(r, t)A,(f, t)dV. (7) 

The 4-vector-potential operator of electromagnetic field is defined by the 
ordinary way as 

t) = £ <fo{Sa*eJ^ < *' + SJ^Se"^. (8) 

a,q 

Here and 6^ - are the operators of creation and destruction of photons 
in states possessed the polarization a(a = 0,1,2,3), momentum q and 
energy u = q. The coupling function g q = ^2tt/u is defined, e.g., in [5]. 
On definition the different unit vectors of polarization e^- are orthogonal to 
each other and obey the normalization conditions (see, e.g., [2]) 

*/i i/ . \ ~> */i i/ ai/ 

9[iv e aq e a l q ~~ 9aa'] 2^, e aq 6 aq9 aa — 9 ■ 

a 

The conditions of the relativistic invariance (see, e.g., ([1],[2],[3])) are sat- 
isfied when the photon operators obey the commutation relations 

Kq, = -gap -A(q-q'); 

Kq, bpq<]- = [V aqt , 6^,]_ = 0. 

Here A(q — q') equal to unity only if q = q' and equal to zero in opposite 
case. It is well known that the non-usual commutation relations for the scalar 
photons lead to the indefinite metric in the photon state space. There are also 
some peculiarities in the definition of physical quantities calculated below. 
By virtue of (4), the equation 

ij- t \t) = H^(t)\t) (9) 

have the exact solution expressed as the direct product of extended coherent 
states [13], namely, 

i*)o = n ex P {-%#) - kqQl q it) + &LA*^)}|°)- ( 10 ) 

a,q 
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Here 

Qaqit) = pqQaq(t), Xaqif) = p\p q Xaqif) (11) 

and ^ 

Qaq(t) = -ig q I dt'e^fM t'y ut \ (12) 
o 

Xaqif) = -- ) {QyflQatf) ~ QUt')Q aq it')}dt'. (13) 



Within our approach an initial vector of state |0) is a direct product of the 
electromagnetic field vacuum state \vac) and the vector \ipo) describing the 
initial state of a particle, |0) = \(fo) (g> \vac) ; below, for simplifications, we 
denote this state also as \ipo,vac) . Assume the vector \ipo) is a free wave 
packet which extends its width in course of time. We must take into account 
that at t > in our interacted system each part of it isn't independent and 
can not be described exactly. E. g., there isn't the wave function for the 
electron but the mostly simplex way to get an explicit information about it 
consists in calculation its density matrix. 

4. The New Perturbation Theory 

If we neglect the corrections producing by the interaction H\ n [ , then the 
equation (10) completely solves the problems with calculating the physical 
quantities of interest. The quantity \Q aq (t)\ 2 has the direct physical mean- 
ing of the mean number of photons possessing the polarization a and the 
momentum q created by the particle to the given instant. Therefore, 

Ak(t) = J2q\Qa q <t)\ 2 (14) 

a,q 

is the mean momentum loss and 

AE(t) = J2^\Q aq it)\ 2 (15) 

a,q 

is the mean energy loss of the particle. Sometimes (e.g., in case of the resting 
particle) the energy loss as well as the momentum loss are zeros. 

In the approximation considered we can calculate also many other physical 
quantities (e.g., the Green function) for the particle with taking into account 
the back influence of electromagnetic quanta on particle's state. The correc- 
tions should be evaluated with the help of a modified perturbation theory. To 
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construct the new perturbation theory with respect to the interaction H-^ t 
we introduce the "zeroth" order evolution operator 



Uo(t) = exp{£Q„j(t)6L- Qls(.t)f>a}-iXas{t)}. 



(16) 



a,q 



The vector (10) now can be rewritten as \t)o = Uo(t)\Q) . Introduce the new 
representation of operators as follows 

A(t) = ui(t)A(t)U (t). (17) 

The vector of state \t) = in this representation obeys the equation 

d 



i-\t) = H^>(t)\t). 



The solution to this equation may be found in a rigorous way via T-exponent 

|0>. (19) 



|i> = Texp{-i/ffS(t')rf*' 





5. Ultraviolet behavior of spectral functions 

The approximate solution to the equation (5) can be found with the help 
of a self-consistent procedure, when we use as \t) an approximate expression 
which functionally depends on f^(q,t) . For the first approximation we sup- 
pose \t) « |t)o • Generally speaking, the quantity f^(q,t) must be chosen in 
a strong connection to the particular problem. We consider now only the evo- 
lution of an one-particle initial state on the background of photon vacuum. 
We begin with evaluation of the left hand of the equation (5). 

Assume ipo is a Gauss wave packet 



(20) 



where A' is a some initial polarization of the electron. Let us represent 
the equation (5) in a more convenient form. Using the well-known Baker- 
Hausdorff formula 

exp(i + B) = exp(i) exp{B) exp(-l/2[i, B]), 

applicable for the case [A, B],A] = [A, B], B] = , we get 



\t) = exp 



a,q \ A j 



10). 



(21) 
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Note, operator Uq is diagonal in spin indices and therefore in our approxi- 
mation only diagonal part of j^(t) contribute to the right hand of (5). The 
straightforward calculation gives: 



(t\jm\t) « (0\Ul(t)j£(t)U (t)\0) « e • exp 



-EIQ^Wf 



a,<f 



X 



a l k-^xMm exp I - E Qa&(*)p&^ 1°)' ( 22 ) 



or, equivalent ly, 



(f « (0|^(^W)|0) = e • exp 
d 3 s 



E \Qa?W 

a,? 



X 



7^ E u Xl p-q/2l^p +m exp 

J{k-p)s 



lie 



p-q/2 



£p+q/2)t] 



X 



E eH*-«-(0| exp - E QmWWz x 



"U-q^k+m exp S QfhMfeibz |o), 

Now, using the relation (see [13]) 

E e ifcf a A £exp | - E QM*)MgW |&, ^ac) = e lk ^\vac p )® 

X,k I ft? J 

expl-EQM^e-^&U l VflC )> 
I ft* J 

where |vac p ) is a particle vacuum state, one can obtain 

f»(q,t) = e ■ E Ux>,?-q~/ 2 Yu\>,p+q/2 exp [i(^-5/2 - £p+q/2)t] x 



(23) 



(24) 



(25) 
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Because Q a q depends on / M , (25) is actually the integral equation for 

Let us assume the initial width of the particle's wave packet is large com- 
pared to the Compton wave length, i.e., 5 ^> 1/mc . In this condition one can 
see in the right hand of (25) the main contribution to the sum occurs from 
the small neighborhood of the point p = ko . Assuming in this neighborhood 
the quantity \p — ko\ is small and comparable with \Ak\ , we can fulfill the 
series development in the power of exponent in the subintegral expression in 
the right hand of (25) with taking into account only terms of the first and 
the second order in (p — ko) . In this case the equation (25) is reduced to 



exp 



Consider the non-relativistic case, when ko <C m , and £^/2+p — £ q/2-p 
£g-/2 • The stationary point of exponent in (26) is found to be 

q t 

$' 2 £q/2 

At sufficiently large width 5 and at q < 2m the last imaginary term in 
the first approximation can be neglected. Moreover, because the width has 
the tendency to increase in course of time (in the non-relativistic case this 
increasing asymptotically develops proportional to time), the last term in the 
right hand of (27) don't increase in course of time. We see the behavior of the 
subintegral function in (26) is defined mostly by the real part of the power 
in exponent. Accepting this assumption, in the first approximation we can 
calculate the sum involved in (26) in the stationary phase approximation and 
get 

t) « e • ux^xfu^+e,* e-^—), (28) 

1 /9 — * — * 

where p m = ko — Ak(ko, t) , e± = (m 2 + (p m ± q/2) 2 ) , Ak(ko, t) - is the 
mean momentum loss of the particle. In this paper we take the equation 
(28) as the basis for constructing the first simple variant of the modified 
perturbation theory. 

6. Mean Electromagnetic Field 

Consider the problem of corrections to the Coulomb field. The mean 
vector-potential of the field is given as A^ 1 = (t\A^\t) . Calculate this vector in 



(26) 



p = k - Ak{t) + i-£ — ; k <C m. (27) 
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the zeroth approximation setting \t) « |0) . With the help of this assumption 
we obtain 

A" = 2e Re{i £ fig^e^-"^ / e^' fM t') p^t') dt'}, (29) 

a,q -oo 

Note, in this expression we have taken into account the negative eigenvalue 
sign for the destruction operator of the scalar photon (see Appendix). Sub- 
stituting (28) and taking into account the summing rules for polarization 
vectors, we get 

^ = 2eRe{iE^™* / u^^u^^x 

a,q -oo 

e i^t' e i(e.-e + )t' p ^ dt 'Y ( 30 ) 

where the function pq{t) is the Fourier component of particle's probability 
distribution in the rest frame. 

— * 

For the uniformly moving particle we have jf = ko and arrive at 



q 

t 



x J e 1 ^' exp{i{e_ - e + )t'} p q it') dt'} : (31) 



— oo 



The potential A^ 1 is represented by the convolution of one for the point par- 
ticle and the probability distribution for the particular state. It is interesting 
to analyze the universal case, corresponding to the point charge, for which 
Pq = 1 . In this case the last formula gives 



q 

or, in 3-representation, 



A» = 2eRe{E g - S ,,,_^,, +g2 eXP 'r;;f;-- )t1 }, (32) 



A = 2e Refe t p wUy/(e+ + m)(g_ + m)+ 
1 q V^+£~ v 



[e+ — m)(e- — m) 2 2/ ■-> e *^* £ 



(/c 2 - g 2 /4 + ^- [fco x g])W ■ ; (33) 

/ UJf? -\- F — F i > 



|^-g/2||A;o + 9/2| ' v + 
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A = 2e Rej^T -j=M== w{, (a(p+ ■ a) + (p_ • a)a)w x > x 

(34) 



exp[igr* — z(£ + — £-)£] 



where p± = (/c ± q/2)y/(e T + m)/{e± + m) . 

Now first we evaluate the scalar potential for the particle found in the rest. 
We substitute ko = and get 

A" = 2e Refe s J»t, _. /2 „ yif/2 ^}. (35 ) 



Because of Uy _ ^/2 u \',q/2 = W £ <f/2 after some algebra the formula (35) 
gives 



2 2mr 

A = e— I K (£)d^ (36) 

where Kq(£) - is the MacDonald's function of zero order. It is important 
to note that at r — > the potential (36) diverges, but only logarithmical as 
(Ziji/tt) ln(e/mr), r <C 1/m. At r > 1/m the potential coincides with 
the Coulomb one. The interaction with the electromagnetic field changes 
particle's form-factor on distances which don't exceed some Compton wave 
lengths. The same behavior manifests the potential which takes into account 
the radiative corrections in the Feinman perturbation theory. 

If we attempt to represent the potential (36) as the usual retarded potential 
corresponded to any charge distribution, then will see, that it is possible only 
if this distribution sufficiently depends on the mass of the particle. The 
dependency of the field on the mass can be connected to the effect of the 
reaction force exited at photon radiation. This circumstance means not the 
charge distribution but the sample action of electromagnetic quanta wrapped 
around the charge creates the field which we observe as the source of the 
Coulomb force. 

Consider the selfenergy for the resting particle. The electromagnetic con- 
tribution to it can be obtained in a usual manner as the energy of an electric 
field 

E = ^j (-VA Q ) 2 dV = ~j ^ A^o dV. 

The integral can be evaluated exactly. First we introduce $ = r Aq . Then, 
using the spherical symmetry and A^Iq = (l/r)d 2 §/dr 2 , bring the energy to 



10 



the form 

8 °r 
E = -r(em) 2 / Kl{2mr) dr = e 2 m. 



In the last expression the value of the integral was used 



oo 2 




/ 



As we see, the electromagnetic contribution to the selfenergy of elementary 
particle is comparatively small. In the case of electron as well as proton the 
electromagnetic correction is equal to the selfenergy multiplied by the fine 
structure constant. We can confirm the famous result the main part of the 
selfenergy for elementary particles has the non-electromagnetic nature. 

Now we should make once more important notation. We can calculate the 
general energy of all the photons wrapping around the particle and found that 
the total its energy exactly consists with the electromagnetic contribution to 
the selfenergy. This result allows us to consider one as a new method to 
calculate the energy of electromagnetic field, which don't use the explicit 
expressions for the electric and magnetic fields. 

— * 

Consider the case of non-relativistic particle. We can set ko « mfo , where 
fo(t) - is the current mean coordinate of the particle, 

mqfo(t) 

c+ — C_ «->-» ; 

£q/2 

e + £- w m 2 (1 + q 2 /4m 2 ). 

In this case we arrive at 

A" = 2e Re{E A m ^^ ~ ^ + ~ e -)t] , (3?) 



and 



A = A s + 2e Re{iJ2^(p-- f^f -)x 

exp[iqf — i(£ + — £-)t] 

UJq + £- — £+ 



(38) 



Note, within the linear in ko approximation the 4-potential (37, 38) obeys 
the Lorentz condition 

dA» _ 
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With the help of simple transformations one can reduce the potential (37) to 
the 3-dimensional integral 

2 \ °r dk °r 
A = -^em / dt' / ^ un / <ircos(/cr)x 

n 2 J J mx (t') I y V y 



l r £fc/2 I 



jif hm\t -t' -T-\f- 



mxo(t' 



£ 



k/2 




K \2m\t -t'-r+ \r - ] \ (39) 

By setting p m = the potential (39) reduces to the potential for resting 
particle (36). 

Within the region r > I /2m the main contribution to the potential (39) 
occurs from the small k < 2m . In this region one can neglect all the nonlinear 
in 7 = k/2m terms in the subintegral function and get 



oo 



A, = - Zmc 2 I dr .y^. T \. \k q (2mc\cr - \r - f (t - r)||) - 
7r J Q \r - r (t - t)\\- 

K (2me(cT+\r-r (t-T)\))\. (40) 

Comparing this result with (39) we see that the previous expression contains 
the additional retardation time r > . The physical explanation of this 
effect can be found on the basis of action of the reaction force appearing at 
the virtual fast k > 2m photon radiation. At this radiation the particle get 
the great (virtual) deflection from the center of the probability distribution. 
This deflections display themselves in (39) in form of the retardation time. 

At the distances r ^> l/2m the MacDonald's functions changes very fast 
compared to the characteristic distance of the potential inhomogeneity. If 
one replaces in (40) the Mac Donald's functions by the corresponded delta- 
functions, the Lienard-Wichert potentials are occur, 



A» = e I dt' - —A 5\t' -t+\r- vt'\] - 8\t' -t-\r- vt'\}}. 
J \r — vt \ 



— oo 



7. Electron Magnetic Momentum Field 

The previous consideration of the mean electromagnetic field in the sec. 6 
has shown, that the Coulomb field isn't the all electromagnetic field of the 
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electron found in the rest. The additional contribution appears from the 
mean value of the vector potential, which general expression is given by the 
formula 

A(r) = 2cRe{E^I w/ 2^/2^^}. (41) 



cj 



After simple algebra one can reduce this expression to 

I(r) = [V$x^], (42) 



where 



^ . N 4m °r sm(x) dx 2 2r f r 
$(r) = — / K = — J K (x) dx , (43) 

7r J Q xVx z + 2mr 7rr J Q 



& el = f^s - is the electron magnetic momentum. Apparently, the expression 
(42) represents himself the vector potential of the electron magnetic momen- 
tum magnetic field. At r — > oo $(r) — ► 1/r, what is right for the point 
magnetic momentum. But at r < m _1 the function $(r) has only the 
logarithmic behavior, simulating the behavior of the main part of potential. 

In the sec. 2 it has been noted, that the interaction with the photon field 
can produce the change in the electron spin direction. Indeed, the explicit 
evaluation of matrix elements (wA',-<f/27 Mw A,<f/2) shows that if in the mostly 
probably processes the condition q ^> ep is fulfilled, then independently of 
any initial polarization in course of evolution the probabilities of all polariza- 
tions become approximately equal. For more explicit proof of this conclusion 
one have to calculate the above matrix elements at ko = : 

t m x 

U \',-q/2 U \q-/2 = ~ "A' A ; 
U \>-q/2® U \,q/2 ' ^q = W \> ° [fy? X ™i W A, (44) 

where - is the unit vector directed along q , and ( j3 = 1, 2 ) - are the 
polarization vectors of transversal photons. Right hand of the second equa- 
tion in (44) occurs in consequence of interaction of the electron magnetic 
momentum with magnetic field of radiating photons. In our approximation, 
when only diagonal part of Hamiltonian included in the ground state defini- 
tion, we don't take into account the spin change during the radiation. In the 
more strict consideration one need to take into account the above speculations 
and to get the regular form of the magnetic momentum field. 
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8. Infrared asymptotics of the photon number 

Consider the asymptotic behavior of the mean number of quanta at 
u = q — > . As it is known, in classical electrodynamics as well as in 
the semiclassical approach there is the characteristic behavior n a q ~ l/u 3 
at u — > . In consequence of this the total number of radiated photons is 
proved to be infinitely large. This result arises due to the rapid increasing the 
photon number at u — > . Now we note that taking into account the reactive 
radiation forces remove the singular behavior of the total photon number. As 
it was shown in [14], the number of radiating photons in collision between 
the moving particle and a resting scatterer is given by the expression 



n a(f (oo) = e 2 gt 



u — qv2 + A u — qv\ + A 



(45) 



where v\ is assuming the velocity before and vi - after collision, \v\ — v<i\ <C 

A = 2e^, fzimif (46 ) 

Here po ~ mv\ . In the non-relativistic limit i;o <C c from (46) we get 
A ~ (4/3)e 2 mt>i . The formula (45) doesn't contain the infrared singularity. 
The deviation from the law n acj "^ 1/u at u — > begins at u ~ A . This 
energy is the less the less the energy of relative motion of the electron and 
the scatterer. 

8. Conclusion 

We see, that problems with divergences in the QED are more the problems 
of the perturbation theory than the QED itself. The divergences character- 
istic for the Feinman-Dyson perturbation theory are eliminated due to the 
physical effects: i) displacement of the particle as the reaction on the own 
radiation, ii) spin overturn processes due to radiation. In the Feinman ap- 
proach the commonly used zeroth order propagators don't take into account 
the electromagnetic interaction at all. Apparently, if we use the more correct 
propagators within the new perturbation theory (which is equivalent to the 
prior "clothing" of particles), the infrared as well as ultraviolet divergences 
at the calculation of physical quantities of interest will be eliminated. 



14 



Appendix 



Coherent States for the Scalar Photon Field 

Because of the non-usual commutation relations for scalar photons, there 
are some peculiarities in formulating the coherent states in this case. Below 
we omit index q in notations of operators <2o,<f for convenience. To construct 
the coherent states for the scalar photon field, applicable within the Lorentz 
gauge and Gupta-Bleuer formalism, we have first to define the Hamiltonian 
for the free scalar photons. In opposite to the transversal and longitudinal 
photons it contains the negative sign in front of the usual expression (see, 
e.g., [3]). 

Ho = -ublbo. 

Consider the Schroedinger equation for the forced oscillator 

ij t \t) = l-ubibo + a(t)bi + a*(t)b }\t). (47) 
The vector of state in the interaction representation obeys the equation 

ij t \t) = [a(t)bl(t)^a^t)b mt). (48) 

Here the interaction representation of operators don't differ from usual, for 
example: bo(t) = b§e~ lut . 

Consider the usual definition of coherent states via D-operator (see, e.g., 

[15]) 

\Q ) = D(Q )\vac) = exp (QA - Q$o)\vac). (49) 
Taking into account bo \ vac) = , we get 

oo f)n 

D(Q )\vac) = e W/2 ■ £ ^(&J)>ac). (50) 

As we see, the sign in the power of the first exponent in the right hand of this 
equation is opposite compared to the ordinary coherent state. The coherent 
state (50) is the eigenvector of the annihilation operator but its eigenvalue 
has the different sign: 

k\Qo) = -QolQo). 

The solution to the equation (48) can be simply obtained with the help of 
famous Feinman-Schwinger approach (see, e.g., in [15]) in the form \t) = 
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e ls °\Qo) • Apply the simple algebra, we get 

• t 

\t) = ex P {^ / Z[Q* (t')Qo(t')}dt'} ■ D(Q (t))\vac), (51) 

to 

where to - is an instant, when we know the initial vacuum state of the photon 
field, 

t 

Q (t) = -iJ a{t')e lLut 'dt'. 

to 
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